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I. Introduction 

Recently Shabat [1] and Spiridonov [2] have discussed potentials which are reflectionless and have 
an infinite number of bound states. In addition, these potentials have the remarkable property of being 
self-similar and can be looked upon as q-deformations of the single soliton solution corresponding to 
the Rosen-Morse potential. Normally, in quantum Lie algebras one takes the underlying space to 
be noncommutative and the deformation parameter q measures deviation from normal analysis. In 
contrast, an interesting feature of Shabat [1] and Spiridonov's work [2] is that they have considered 
the same problem in a commutative space with the q-deformation arising from the specific nature 
of the potential. Spiridonov [3] has also considered the deformation of parasupersymmetric quantum 
mechanics [4] and obtained potentials which can be regarded as a q-dcformation of the two soliton 
solution corresponding to the Rosen-Morse potential. 

Another interesting advance of recent years has been in supersymmetric quantum mechanics [5], 
where new insight into exactly solvable potentials has been obtained through the concept of shape 
invariance. It has been shown [6] that for supersymmetric partner potentials V±(x,a ) satisfying the 
properties of shape invariance and unbroken supersymmetry, one can write down the energy eigenvalues 
algebraically. Subsequently it was shown that both the eigenfunctions [7] and the scattering matrix [8] 
can also be obtained algebraically for these potentials. The shape invariance condition is given by 

V + (x,a ) = V-(x,a 1 ) + R(a ), (1.1) 

where ao is a set of parameters, a\ = /(ao) is an arbitrary function describing the change of parameters 
and the remainder i?(ao) is independent of x. Certain solutions to the shape invariance condition are 
known [9] including essentially all the standard problems discussed in quantum mechanics textbooks. 
In all these cases a\ and ao have been related by a translation (oi = ao + a). Careful analyses with this 
ansatz have failed to yield any additional shape invariant potentials [10]. Indeed it has been suggested 
[11] that there arc no other shape invariant potentials. Although a rigorous proof has never been given, 
no counter examples have so far been found either. 

In this paper we show that the Shabat-Spiridonov (SS) self-similar potentials can be understood 
within the framework of shape invariance. In particular, we show that by using a scaling ansatz 
(ai = qao) for the change of parameter ao, a large class of new shape invariant potentials which are 
reflectionless and possess an infinite number of bound states can be obtained [12].* Our potentials 
contain the self-similar potentials of Shabat [1] and Spiridonov [2], but are considerably more general. 

* This is slightly misleading in that a reparameterisation of the form a\ = qao can be recast as 
o'l = a 'o + ol merely by taking logarithms. However, since the choice of parameter is usually an 
integral part of constructing a shape invariant potential, it is in practice part of the ansatz. For 
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The plan of the paper is the following. In Section II, we discuss in some detail the self-similar 
potentials of Shabat and Spiridonov (SS). An unfortunate feature of these potentials is that they are 
not known in analytical form for all x and we therefore give graphs of these potentials for some values 
of the deformation parameter, which we denote by p (0 < p < 1). 

In Section III, we briefly discuss the shape invariance condition within the formalism of supersym- 
metric quantum mechanics. Using a scaling ansatz (a\ = qa ), we obtain a large class of new reflection- 
less shape invariant potentials. Explicit expressions for the eigenvalues, eigenf unctions and transmission 
coefficients for these potentials are derived. These potentials can be viewed as q-deformations of a one 
dimensional harmonic oscillator or of the single soliton solution corresponding to the Rosen-Morse 
potential. The self-similar potentials of SS [1,2] are rederived as a special case. 

In Section IV, we discuss a new technique that essentially solves the inverse scattering problem 
for this wider class numerically and so enables us to calculate these potentials. Examples of the results 
obtained are discussed. 

In Section V, we give the Taylor series expansion of the potentials for large x. By using these 
reflectionless potentials as solutions of the KdV equation, we then estimate the area under them and 
indicate how one can also estimate higher moments for these cases rigorously even though the potentials 
are not known in analytical form. Using the ground state wave function for these potentials, we also 
give graphs of a continuous parameter family of potentials which are strictly isospcctral to one of the 
self-similar potentials. 

Section VI introduces multi-step scaling ansatze for the change of parameters and hence obtains 
new shape invariant potentials which can be looked upon as q-deformations of the multi-soliton solu- 
tions corresponding to the Rosen-Morse potential. 

In Section VII, we discuss various other ansatze for connection between parameters a\ and an 
and obtain yet more new shape invariant potentials. Explicit expressions for the eigenvalues and the 
eigenfunctions of these potentials are also given. 

Finally, in Section VIII, we summarize the results of this paper and indicate some open problems. 

II. Self-Similar Potentials 

Shabat [1] and Spiridonov [2] discussed an infinite chain of reflectionless Hamiltonians given by 

instance, in section III, we will construct potentials by expanding in do, a procedure whose legitimacy 
and outcome are clearly dependent on our choice of parameter (and hence reparameterisation) . Note 
that, although the construction is non- invariant, the resulting potentials will still be invariant under 
redefinitions of a . 
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(h = 2m= 1) 

H n =P 2 + V n (x), n = 0,1,2,... (2.1) 

with 

V (x) = W 2 -W^ + Co (2.2) 

and 

V n+1 {x) = V n (x) + 2W' n (x). (2.3) 
The various superpotentials W n {x) satisfy the following set of differential equations 

Wl + W' n = W 2 n+l - < +1 + C n+1 , n = 0, 1, 2, . . . (2.4) 

where the C„ are arbitrary positive constants. It is amusing to note that relations (2.4) arise naturally 
in the framework of parasupersymmetric quantum mechanics [4] . Let us assume at this stage that all 
W n (x) are such that the functions 

4 n) (x)c<exp[-y W n (y)dy], (2.5) 

are square integrable, and hence correspond to the ground state wave function of the Hamiltonian H n . 
In this case one has the standard situation of unbroken supersymmetry and the potential Vi + \{x) has 
one bound state less than Vi(x). Using eqs.(2.1) to (2.4), it follows that the eigenvalues of H are given 

by 

m 

^ 0) =E C - m = 0,1,2,... (2.6) 

i=0 

while the corresponding eigenfunctions are given by 

^(x)a(P + iW Q ){P + iWi) (P + iW m ^)^ n \x). (2.7) 

It should be noted here that the Hamiltonian Hj has the same spectrum as H except that the lowest 
j levels of Ho are missing. 

In general, it is not possible to determine these potentials unless one imposes some extra con- 
straints. SS specify superpotentials by demanding that all superpotentials W n {x) satisfy the self-similar 
ansatz 

W l {x)^p l W{p l x) (2.8) 

with < p < 1. Thus, there is just one unknown function W(x). On using i = and 1 in eq.(2.8) and 
eq.(2.4), one obtains the following finite-difference differential equation defining W(x): 

W 2 (x) + W'(x) = p 2 W 2 (px) - p 2 W'(px) + d (2.9) 
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where the prime denotes differentiation with respect to the argument. Eqs.(2.8) and (2.9) are the key 
statements underlying the concept of self-similarity. On using i = 2, 3, . . . and eqs.(2.4), (2.8) and (2.9) 
one then concludes that 

C n = (p 2 )"- 1 ^ (2.10) 



and hence the mth eigenvalue of H is given by 

(1-P 2 ) 



f 71 _ n 2ra) 

EM=Ca+ ^ * , m = 0,l,2,... (2.11) 



One can choose the arbitrary constant Co to be zero, which corresponds to taking E^ = 0. An 
alternate convenient choice is to pick Co such that 

Lim^ooE^ = 0, (2.12) 

which gives 

£ ™=-(rr7)' m = °M--- ( 2 - 13 ) 

One can try to find W{x) by solving the finite-difference differential eq.(2.9) in a Taylor series 
form near x = 0; if 

oo 

W(x) = b 3 xJ ( 2 - 14 ) 

3=0 

then 

b 1 (l+ P 2 ) + bl(l-p 2 )=C 1 (2.15) 

and 

b ^ = - (1 +i)(C%) i bmbj - m ' 3 = i; 2i • • • (2 - i6) 

Normalizability of wave functions is ensured if VF(x) is a continuous function, positive at x — > oo and 
negative at x — > — oo. This is the case if one chooses &o = 0. In this case, it follows from (2.16) that 
all even coefficients bj (j = 0, 2, 4, . . .) are zero and hence W(x) = —W(—x). In particular one finds 

Some special cases are worth noting. At p = 1 the solution of eq.(2.9) is the standard one- 
dimensional harmonic oscillator with W(x) = C\x/2, while in the limit p — > 0, the solution of (2.9) is 
the one soliton superpotential corresponding to the Rosen-Morse potential given by 



W(x) = V / Citanh( v /Cix). 
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(2.18) 



Hence the general solution to eq.(2.9) with < p < 1 can be looked upon as a deformation of the 
hyperbolic tangent function with p acting as the deformation parameter. Notice that the number of 
bound states as given by (2.13) increases discontinuously from just one at p = to infinity for p > 0. 

It is important to note that the superpotentials W(x, C\) which solve the self-similarity condition 
(2.9) have the simple scaling property 

W{x,Ci) = y/ClF{y/Oix). (2.19) 

Thus, for a particular p one only needs to find W(x, C\) for any one non-trivial value of C\. Knowing 
W(x), one immediately knows V (x) and the other potentials V n {x) can be recursively obtained from 
Vq(x) using eqs.(2.1) to (2.4), (2.8) and (2.10), so that the whole chain of potentials are known in 
principle once Vq{x) is known. 

Unfortunately, merely knowing the Taylor series about x — is not sufficient if one wants to carry 
out this programme, since simple arguments show this series (2.17) to have a radius of convergence R, 
where 

l<VCiR<\<j)±^, 0<p<l. (2.20) 

The lower bound derives from noticing that the series coefficients are smaller for < p < 1 than 
for p = 0, the latter being essentially the expansion of tm\\\{^/C\x) and known to have radius of 
convergence 7r/2 (see also [13]). The upper bound involves realising that there still has to be a pole on 
the imaginary axis when p > 0: from (2.17) one sees that along that axis W(ix) = iuo{x) inside any 
radius of convergence, where u(x) is a real function satisfying 

-L0 2 (x) + LU'(X) = -p 2 L0 2 ( P x) - p 2 L0'( P x) + C U (2.21) 

but also that lo(x) > u)(px) and ui'(x) > u'(px), so that 

lo\x)(1 + p 2 ) > lu 2 (x)(1 - p 2 ) + & (2.22) 

implying that ui(x) grows faster than tan(v / C i i\/l — p 2 x/^/l + p 2 ) and so has the requisite singularity. 

In the absence of a solution to (2.9) in terms of elementary functions for all p, one must resort 
to some sort of numerical determination other than by trying to sum the Taylor series. The most 
direct approach relies on noticing that since px < x, if one already knows W(x) in an interval, then 
one knows the right hand side of (2.9) in a larger interval. This allows us to treat the right hand side 
as an already known function of x and thus integrate up the equation using a (fourth order) Runge- 
Kutta method. To start this process off, the Taylor series was summed to 75 terms on an interval 
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well inside the radius of convergence. The superpotentials thus obtained have been checked both by 
comparing them to the summed series throughout the region where the latter is still valid and by 
direct numerical integration of the Schrodinger equation with the spectra checked against (2.13). In 
both cases the degree of agreement is extremely high. Examples of the kind of superpotentials and 
potentials obtained are shown in Figures 1 and 2. Note that, having independently calculated these 
functions, we see no evidence of the oscillations in them reported by in [13]. 

Aside from its practical utility, the insight underlying the numerical determination also strength- 
ens confidence that (2.9) actually has a solution. Because any solution in a finite interval can be 
analytically continued to arbitrarily large x, the question of existence clearly reduces to establishing it 
in a neighbourhood around x = 0. But this is precisely the place where a convergent Taylor series is 
known to exist and this is sufficient. 

Additional analytic properties of these potentials will be described in Section V, but first we wish 
to introduce our more general class of potentials to which the methods of that section will also apply. 

III. Shape Invariance With Scaling Ansatz 

A fresh impetus to the study of exactly solvable problems in nonrelativistic quantum mechanics 
was provided by Gcndcnshtein [6] with the introduction of shape invariant partner potentials within 
the framework of supersymmetric quantum mechanics. To set our notations, we give a quick review 
of both supersymmetric quantum mechanics and shape invariance [9] . The partner Hamiltonians H± 
are given by 

H- = A + A, H+ = AA + , (3.1) 

where (h = 2m = 1) 

A=4-+W(x), A+ = 4--W(x), (3.2) 
ax ax 

so that the two partner potentials V±(x) can be expressed in terms of the superpotential W(x) thus 

V±{x)=W 2 {x)±W'(x). (3.3) 

From here it follows that all the energy eigenvalues of H± are positive semidefinite. Further, it turns 
out that in case SUSY is unbroken the ground state energy of one of the two Hamiltonians is zero and 
all other energy eigenvalues of H± are paired 

4 -) =0, =£<+>. (3.4) 

Here, for convention's sake, we always consider the situation of unbroken SUSY and so the ground 
state energy of H- is zero. The corresponding eigenfunction ipQ~\x) (which satisfies Atj)Q~\x) = 0) 
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turns out to be 

^\x) = Ne-f W{y)dy . (3.5) 

One can also show that because of SUSY the eigenfunctions and scattering amplitudes of the two 
partner Hamiltonians are also related 

= A^Jy/s^, - A + ^/4W\ (3.6) 

W 4- ik 

R - {k) = mr^ik R+{k) > (3 - 7) 

T -^ = W7^k T+(fc) ' (3 - 8) 

where k = (E — Wl) 1 ' 2 and k' = (E — W 2 ) 1 ' 2 with W± = W(x = ±oo). 

If the pair of SUSY partner potentials V±(x) defined by eq.(3.3) differ only via the parameters 
that appear in them, then they arc said to be shape invariant [6]; that is, if the partner potentials 
V±(x, ao) satisfy the condition (1.1). In terms of the superpotential W, this shape invariance condition 
reads 

W 2 (x, a ) + W'(x, ao) = W 2 {x, oi) - W( x,ai) + R(ao). (3-9) 

The common x-dependence in V- and V+ allows a full determination of energy eigenvalues [6], eigen- 
functions [7] and scattering matrices [8] algebraically. One finds 

n-l 

E { -\a ) - Y, R ^ ak ^ E o~ ] M = 0, (3-10) 

fc=0 

tp^~\x,a ) = A + (x,a )A + (x 1 ai) . . . A + (x, a„_i)?/; _) (x, a n ). (3.11) 

It is still a challenging open problem to identify and classify all the solutions to the shape invariance 
condition (3.9). Certain solutions to it are known [9] and they include essentially all exactly solvable 
problems discussed in standard texts on quantum mechanics. For all these, a\ and ao are related by 
a translation. Careful analysis with this ansatz has failed to uncover any additional shape invariant 
potentials [10] and in fact it has been suggested that there are no others [11]. We shall now show 
that this is not the case since a large number of new shape invariant potentials can result from a new 
scaling ansatz 

ai = ga (3-12) 

where < q < 1, a choice motivated by recent interest in q-deformed Lie algebras. Our approach 
includes the self-similar potentials discussed in the previous section as a special case. 



Consider an expansion of the superpotential of the form 

oo 

W(x,a ) = ^2gj{x)a 3 . (3.13) 

3=0 

Using eqs.(3.12) and (3.13) in the shape invariance condition (3.9), writing R(a ) in the form 

oo 

R(ao) = Y,Rj< ( 3 - 14 ) 

3=0 

and equating powers of do yields 

2g' (x) = Rq, g[(x) + 2d 1 g (x)g 1 (x) = ndi, (3.15) 

n-l 

g' n (x) + 2d n g {x)g n (x) = r n d n - d n g j {x)g n - j {x) (3.16) 

i=i 

where 

r n = R n /(l-q n ), d n = (l-q n )/(l + q n ), n = 1,2,3,... (3.17) 

This set of linear differential equations is easily solvable in succession to give a general solution of 
eq.(3.9). Let us first consider the special case 

9o(x) = (3.18) 
which implies i?o = 0. The general solution of (3.16) is then 

/n— 1 
dx r n -^2g j (x)g n - j (x) , rc = 2,3, ... (3.19) 
i=i 

where without any loss of generality we have assumed all the constants of integration to be zero. The 
shape invariance condition thus essentially fixes the g n (x) (and hence W(x,ao) via (3.13)) once R(ag) 
is specified, i.e. once the set of r n are chosen. Implicit constraints on this choice are that the resulting 
ground state wavefunction (3.5) be normalisable and, so that the spectrum (3.10) is sensibly ordered, 
that R(q n a ) > 0. 

The simplest case is r\ > (positivity required to ensure normalisable wavefunctions) and r n = 0, 
n > 2. Here (3.19) takes on a particularly simple form g n {x) = PnX 2 "^ 1 with 

d ri—l 

/?i= din, In = - (2w " t) S ftA-j" = 2, 3, . . . (3.20) 

and so 

oo 

Wfo a ) = Aoo* 2 *" 1 = V^^(V^^)- (3-21) 
9 



For a\ — qao this now gives 

W(x, 01) - VqW(y/qx, oo), (3.22) 
hence in this special case the shape invariance condition (3.9) becomes 

W 2 (x, oo) + W'fo oo) = gW*(v^s, oo) - ^^^V^, «o) + aori(l - q). (3.23) 

Comparing this to (2.9), one thus sees that the case r n = 0, n > 2 corresponds to the self-similar W 
of Shabat and Spiridonov provided one writes j 2 = diria and q = p 2 . In fact, if instead of choosing 
r n = 0,n > 2, any one r„ (say r 3 ) is taken to be nonzero and qi is replaced by p 2 then one again 
obtains the self-similar potentials. In these instances the results obtained from shape invariance and 
self-similarity are entirely equivalent and the Shabat-Spiridonov self-similarity condition turns out to 
be a special case of the shape invariance one (3.9). 

However, it is necessary to emphasize here that shape invariance is a much more general concept 
than self-similarity. For example, if we choose more than one r n to be nonzero, then shape invariant 
potentials are obtained which are not self-similar. As an illustration, consider r n = 0, n > 3. Using 
eq.(3.16) one can readily calculate all the g n (x), of which the first three are 

gi(x) = dxnx, g 2 (x) = d 2 r 2 x - ^dfr^x 3 , 

2 2 
m{x) = --d 1 r 1 d 2 r 2 d 3 x 3 + —d\rld 2 d 3 x h . (3.24) 
6 15 

Notice that in this case W{x) contains only odd powers of x. This makes the potentials V±(x) 
symmetric in x and also guarantees unbroken SUSY. It is convenient to define the combinations 
Ti = G?iria = 7 2 and Y 2 = d 2 r 2 a^. Then, the energy eigenvalues follow immediately from eq.(3.10) 
and (3.14) (0 < q < 1): 

F (-) (n) - r (Ufa!) , r (i + g 2 )(i-g 2 ") „ _mo r o 2 ^ 

K'\ a o)-Li (T~qj 2 (l - q 2) ' n-0,1,2,... (3.25) 

while the (unnormalized) ground state wave function is 



x, ao) = exp 



x 2 



-:r(ri + r 2 ) + — (d 2 r( + 2d 3 r 1 v 2 + d 4 r 2 ) + o(x b ) 



(3.26) 



The excited wave functions can be recursively calculated using eq.(3.11), though usually it is more 
convenient to use the relation 

^-\x,a ) = A + {x,ao)^-\{x, ai ). (3.27) 
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We can also calculate the transmission coefficient of this symmetric potential (k = k') by using 
relation (3.8) and the fact that for this shape invariant potential [8] 

T + (k, ao) = T-(k, oi = qa ). (3.28) 

Repeated application of eqs.(3.8) and (3.28) gives 

[ik - W(oo, a )}[ik - W(oo, ai )}....[ik - W(oo, a ra _i)] 
1 ' oJ ~ [ik + W(oo, a )][ik + W(oo, oi)]....[ifc + T^(oo, o„_i)] 1 ' " j 1 ' 

where 



VF(oo, Oj) = ^EU~Ep. (3.30) 

Now, as n — > oo, a„ = a"&o — > (0 < a < 1) and, since we have taken go(x) = 0, one gets 
W(x; a n ) — ► 0. This corresponds to a free particle, for which the transmission coefficient is unity; as a 
result the reflection coefficient R-(x; ao), vanishes and the transmission coefficient is given by 



T_(fc,a ) = n 

3=0 



ik — W(oo, a,j) 



ik + W(co, dj) 



(3.31) 



Note that this is a general expression for T derived using only a„ — q n a and the fact that V-(x) is a 
symmetric potential (k = k'). In particular, no assumption has been made regarding the values of r n . 
In the case when only one of the r„ (say n) is non-zero, eq. (3.31) gives the transmission coefficient 
of the self-similar potentials (3.13) to (3.16). 

The above discussion keeping only n , r 2 ^ can readily be generalized to an arbitrary number of 
nonzero rj. The energy eigenvalues for this case are given by (Tj = djTja 3 ) 

^M-E ^ 1 ^^^ - , = 0,1,2,... (3.32) 

3 

All these potentials are also symmetric and reflectionless with T_ given by eqs.(3.31). 

The limit q — > is particularly simple and again yields the one-soliton Rosen-Morse potential with 
W = at&nhax. Thus results corresponding to different choices of R n can be regarded as multiparam- 
eter deformations of this potential. 

Finally, let us consider the solution to the shape invariance condition (3.9) in the case where i?o 
is nonzero, so that go(x) = ^RqX from (3.15) rather than being zero. One can again solve the set of 
linear differential eqs.(3.16) in succession using this go(x) and hence obtain gi(x), g2(x), ■ ■ ■ Further, 
the spectrum can be immediately found using eqs.(3.10) and (3.14); for example, in the case of an 
arbitrary number of nonzero Rj (in addition to Rq), the spectrum is given by 

^n^r /'+ff-^' (3.33) 

3 
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which is the spectrum of a q-deformed harmonic oscillator [14]. It should be noted here that, unlike 
the usual q-oscillator where the space is noncommutative, but the potential is normal (w 2 x 2 ), in our 
approach the space is commutative, but the potential is deformed, giving rise to such a multi-parameter 
deformed oscillator spectrum. 

IV. Numerical Results 

Explicit determination of the SS potentials (described in Section II) crucially depended on the 
scaling property 

W(x, a ) = yf^F(^x), (4.1) 

displayed by the solutions, which allowed W'(x,a n ) to be related to W(y/qx, ao) instead of merely 
W(x, ai) = W(x,qao). However, such scaling is not a property of the solutions in Section III when 
more than one r n is non-zero, as can be seen from the series expansion (3.24). When only r\ and r 2 
are non-zero, there is a generalisation of the form W(x, ao) — y/riaoF(^naox, r20o/ri) which relates 
the behaviour at x to that of another problem (that corresponding to calculating W(x,a ) with a 
different r 2 ) at y/qx. By forming a ladder of these potentials (in which r 2 is tending to zero and hence 
the problem towards the special SS case), it should be possible to determine W(x, ao) using this fact. 
However, we chose to devote this section to a method that emphasizes W(x, a) as a function of both 
x and a and which more readily generalises to arbitrary r n . 

Intuitively one would still expect the series (3.13) to be convergent for either iora sufficiently 
small. Thus in the (x, a) plane we can assume that W(x, a) is calculable to arbitrary accuracy close 
to either axis and the problem reduces to continuing this knowledge out into the plane. Defining sum 
and difference functions 

S(x, a) ee W(x, a) + W{x, qa) (4.2a) 
D(x, a) ee W(x, a) — W(x, go), (4.26) 

for ai = gao, the shape invariance condition (3.10) becomes 

^ = -S(x,a)D(x,a) + R(a) (4.3) 
ax 

where R(a) is a known function. Now if one knows W(x, a) for x < X, a < A, an Euler step (since D is 
only known in x < X , there's not enough information available for that to be a Runge-Kutta one) will 
give S(X + h,a), again for a < A. One must now invert (4.2a) to convert this knowledge of S(X + h, a) 
into information about W(X + h, a). Iterated use of (4.2a) relates W(X + h, a) to W(X + h, q n a) and 
n known values of S. For some sufficiently large n, W(X + h, q n a) can be calculated using the Taylor 
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series (3.14) and so one can indeed determine W(X + h, a) for all a < A. Breaking the (x, a) plane up 
into a grid and using some suitable interpolation method for the points in between, one can iterate this 
Euler step up through x and so obtain W(x, a) numerically for values of x and a limited by computing 
constraints only. The only inputs are R(a) and the series approximations for small x and a. 

A program to implement this scheme has been developed and its results for the special SS case 
were shown to agree very well with the earlier (more accurate, but overly specialised) program. As 
an example of the new potentials this permits us to consider the potential corresponding to n = 1, 
r 2 = —1 with all other r n = 0. Provided a < 1/(1 + q) the spectrum given by (3.25) is well-ordered. In 
Figure 3 we display the potential calculated for q = 0.3 and a = 0.75, along with its partner potential 
and the exact spectrum found from (3.25). (These eigenvalues have also been checked numerically). 
Note that for this choice of parameters, V- (x) is a double well potential, whereas its shape invariant 
partner V+(x) is a single well. This situation is unlike the (non-shape- invariant) examples discussed 
in Ref. [15], where the SUSY partner of the initial double- well potential has a sharp 5-like spike at its 
center. Apart from being the first shape invariant double-well, this example stretches naive intuition 
concerning shape invariance. Furthermore, this example barely indicates the variety of behaviour 
available by altering q, a and the r„ in this new class of potentials. 

Finally, we note that the basic idea of this section can be divorced from the details of the Taylor 
series. The restriction to symmetric potentials gives W(0, a) = 0, to be used as a boundary condition 
for the intial Euler step. To invert (4.2a) one can use the infinite series 

W(x, a) = S(x, a) - S(x, qa) + S(x, q 2 a) - S(x, q 3 a) + ... (4.4) 

which is convergent provided W(x, a) — > as a — > 0. More useful numerically is the fact that, because 
it is alternating, this series can be truncated with a rigorous bound on the error and without needing 
to calculate some W(X + h, q n a) by use of a Taylor series. Neither is the reliance on oi = qa a terribly 
restrictive: one can always redefine the parameters to obtain this. The only crucial constraint is 
that the method applies exclusively to symmetric potentials holding infinitely many bound states and 
corresponding to a chain of Hamiltonians (H-, H + etc.) which tends asymptotically towards a free- 
particle one (i.e. W(x, a n ) — » as n — > oo). Otherwise the sole input is R(a) (expressed in terms of the 
appropriately defined a), which in general should be deducible from any desired (possible) spectrum. 

V. Analytic Results 

Although undeniably useful, simply being able to calculate W is an unsatisfactory state of affairs 
unless there is also a body of complementary analytic results. This section therefore brings together 
several approaches by which such results can be gathered. 
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In our method of constructing the potentials, the Taylor series about x = played an essential 
role (as it did in [1] and [2]). However, to get a better insight into the potentials, it may be worthwhile 
to also know the Taylor series of W(x) around x — > oo. For simplicity, we now restrict our attention 
to the SS family and return to using the original parameter p = ^fq. Substituting t = 1/x in eq.(2.9) 
yields 

W 2 (t) - t 2 W'(t) = p 2 W 2 (t/ P ) + t 2 W'(t/p) + d. (5.1) 



On substituting* 



in this equation one obtains 



and 



3=0 



a 2 = C'i/{l-p z ), oi =0 (5.3) 



a, =-(j- 1) (^U) ^ - ^ E a ™ a — i = 3, 4, . . . (5.4) 
Thus we find that asi^ +oo 

^)^/tS + V-(^H + --- (^) 
y(l-p 2 ) aoVl-p/a;' 3 

where a 2 is an arbitrary constant. This arbitrariness is due to the fact that W(0) — has not been 
imposed while deriving (5.5). In fact it is not easy to do so since the series (5.5) is valid for large x. 

It has already been established that we are dealing with reflectionless, symmetric potentials for 
which the infinite spectra of eigenvalues arc known exactly in closed form. This type of problem has 
already been well studied, but the standard inverse scattering method [16] has proved too cumbersome 
to be of much practical use in deriving these potentials in this non-trivial context. However, certain 
well-known, related results can be used to quite strongly constrain the potentials: the point is, being 
reflectionless, these can be regarded as a solution of the KdV equation at time t = [17]. Now it is 
known that such a solution as t — > ±oo will break up into an infinite number of solitons of the form 
2k 2 sech 2 kiX. On using the fact that KdV solitons obey an infinite number of conservation laws corre- 
sponding to mass, momentum, energy etc., one can immediately obtain constraints on the reflectionless 
potentials by using the known solutions at t — ► ±oo. For example, the first three conservation laws are 

/ V (x)dx = Y^ / V®(x)dx (5.6) 



* We assume, in the absence of evidence to the contrary, that the simple Taylor series is the appro- 
priate expansion. 
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V 2 (x)dx 



OO 

OO „ 

£ / [VW(x)] 2 (to 



l,dK«(a;) )2 - 



da; 



where 



V l {x) = -2kfsech 2 kiX 



(5.7) 



(5.8) 



(5.9) 



with hi = k a p l and A; 2 = Ci/(1 — p 2 )- Using eq.(5.9) it is straightforward to evaluate the right hand 
sides of eqs.(5.6) to (5.8) and we find 



oo ^ 

/ Vo(x)dx = -4j2h = ^ 



(5.10) 



oc 

/ 



i=0 



3 (l-p3) 



*>4(f) 2 ] 



32 v ^ , c 32 /en 



5 (1-p 5 ) 



(5.11) 



(5.12) 



thereby providing strong constraints on the potential Vo(x). 

All deformations of potentials considered so far have been such that the spectra obtained are q (or 
p) dependcndcnt. Before ending this section, it is worth remarking that, as with any potential, there 
are also distortions of the V n (x), with deformation parameter A, which leave the spectra unchanged. 
Using the techniques of supersymmetric quantum mechanics, one can construct a large class of strictly 
isospectral potentials. For example, using any one of the Wi(x) as given by eqs.(2.5), (2.8) and (2.14) to 
(2.17) one can immediately obtain a one parameter family of strictly isospectral reflectionless potentials 
V n (x, A) by using the formula [8] 



V n (x, A) = V n (x) - 2— ln(J n (a;) + A) 



(5.13) 



where V n (x) can easily be obtained using eqs.(2.1) to (2.3), (2.8) and (2.14) to (2.17), A is any arbitrary 
parameter (A > or A < —1) and 

'[4 n) (y)} 2 dy. 



In(x) 



(5.14) 



Here is as given by eq.(2.5) which can be explicitly obtained by using eqs.(2.8) and (2.14) to 
(2.17). As an illustration, we give graphs of the Vq(x, A) obtained from the SS potential with p = 0.5 
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for various values of A in Figure 4. Large values of A correspond to the original SS potential (see Fig. 4c). 
As A takes on values closer to zero, the potential gradually breaks into two pieces, one corresponding 
to the E — state only and the other containing the remaining energy levels [8] . 



VI. Shape Invariance in More Than One Step 

Having obtained potentials which are multiparameter deformations of the one soliton solution of 
the Rosen-Morse potential, an obvious question to ask is if one can also obtain deformations of the 
multi-soliton solutions. The answer is yes and as an illustration we now explicitly obtain multiparam- 
eter deformations of the two soliton case. The desired deformation is achieved by extending the usual 
shape invariance ideas to the more general concept of shape invariance in two steps. 

Consider the unbroken SUSY case of two superpotentials Wo(x,a ) and Wi(a;, flo) such that 
are the same up to an additive constant. 

V {+ \x 7 a ) = V 1 { -\x,a ) + R(a a ), (6.1) 

or equivalently, 

W$(x, ao) + W^x, ao) = W?(x, ao) - W[(x, a ) + R{a ). (6.2) 
Shape invariance in two steps means that 

v} + \x,a ) = Vt\x iai ) + R(a ), (6.3) 

that is 

W?(x, a ) + W[(x, a ) = W^x, oi) - Wfa, oi) + R(a ). (6.4) 

For the above situation, the energy eigenvalues and eigenfunctions of the potential Vq~\x, a ) can be 
algebraically calculated as shown below. 

Unbroken SUSY implies zero energy ground states for the potentials Vo~^(a;, do) and v} \x, ao): 

E { Q - )a = 0, E ( - ]1 = 0. (6.5) 

The degeneracy of energy levels for supersymmetric partner potentials yields 

4 +)0 M = Ei-\°(a ), E^(ao) = E { ~^(a ). (6.6) 

From cq.(6.1) it follows that 

Ei^(ao) = E^(ao) + R(ao). (6.7) 
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For the special case n = 0, eqs.(6.6) and (6.7) give 

E ( - )0 = R(a ). (6.8) 
Also, the shape invariance constraint (6.3) gives 

E^ 1 (a ) = E^°(a 1 ) + R(a ). (6.9) 
Using eqs.(6.6), (6.7), (6.9) and some algebra, one gets 

E { n+i M = (oi) + ^(ao) + £(«(>)■ (6-10) 
These equations can be solved recursively to get 

n-l 



4 )0 = E[%)+%)1. 

fc=0 
n-l 

E { 2 n+i = E [^(Ofc) + R{a k )] + R(a n ). (6.11) 

The above discussion has been completely general and is valid for any change of parameters, a\ — f(ao). 
Following the treatment of Section III, we now take the scaling ansatz a\ = qa and expand the 
superpotentials Wo and W\ in powers of a . 

oc 

W„(x,a ) = J2 9j(x)4 (6.12) 
3=0 

OO 

W 1 {x,a )=J2 hj(x)a&. (6.13) 

3=0 

Further, write R and R in the form 

R{a ) = jr Rja , R(a ) = Rja{. (6.14) 

3=0 3=0 

Using these in eqs.(6.2) and (6.4) and equating powers of a yields (n = 0,1,2,. . .) 

n n 

9n + 9j9n-j = h 3 h n-j ~ h 'n + Rn (6.15) 

3=0 j=0 

n n 

K + h 3 h n-j =q n J2 9j9n-j ~ q n 9'n + Rn- (6.16) 
3=0 j=0 
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This set of linear differential equations is easily solvable in succession. Let us first discuss the special 
case 

go(x) = h (x) = 0, (6.17) 

which implies that Ro = Ro = 0, and further assume that R n — R n — 0,n > 3. In this case one can 
readily calculate all g n {x) and h n (x); the first two of each are 



9i = 



( Ri ~ Ri ) 
(1-9) 



R 2 x 



-x, 92 = 



Cgg - R 2 ) 



x + 



3(1 - g) 3 



(l-q){Rl-R\)-2{l + q)R 1 R 1 



(R 1 -qR 1 ) 

hi = —r r—X, 

(1-9) 

(1 + q)R\ + (1 + q)(l - q 2 )R\ - 2q{\ - q)R 1 R 1 



m + q)(l-qf^'^'^'^ — — J" (6J8) 

ft may be noted that both Wq and W\ contain only odd powers of x so that the potentials and 
arc all symmetric in x and SUSY is unbroken. The energy eigenvalues can be obtained from 
eqs.(6.8) and (6.11). 

(6.19) 
(6.20) 



E[ )0 (a )=Ria + R 2 al 

2 

eL )0 M = J2(Rj + i^ (^O, 



and 



(6.21) 

V I — fl> 

3 = 1 " J=l 

For the special case when R 2 = R2 = 0, the spectrum has been obtained previously by Spiridonov 
from consideration of self-similar potentials [3]. However, the spectrum in the general case given by 
eqs.(6.20) and (6.21) cannot be obtained in such a fashion. The energy eigenvalues E^ 1 are now 
immediately obtained from eq.(6.9) and the energy cigenfunctions and transmission coefficient for 
these rcficctionless potentials can also be found using eqs.(3.5), (3.27), (3.30) and (3.31). Further, the 
above discussion can be readily generalized to an arbitrary number of nonzero Rj , Rj . 

The limit q — > of the above equations is particularly simple and yields the two soliton solution 
of the Rosen-Morse potential, i.e. 



W = 2 tanh \f~kx, W x = V^tanh VlL 



(6.22) 



provided R — 3R. Thus our results can be regarded as multi-parameter deformations of this potential. 

Finally, it is clear that one can easily generalize this procedure and consider shape invariance with 
a scaling ansatz in 3,4,. . .p steps and thereby obtain multi-parameter deformations of the 3,4,. . .p 
soliton Rosen-Morse solution. 
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VII. Shape Invariance with a Non-Scaling Change of Parameters 

We have so far obtained new shape invariant potentials for ax and a related by the scaling ansatz 
(ai = qa ). Are there shape invariant potentials where ax and a are neither related by scaling nor by 
translation (ax = a + a)l We now demonstrate the existence of yet other possibilities by obtaining 
potentials for ax = qa^ and ax — qao/{l + pao)- 

First consider the case when 

ax = qa p (7.1) 

where p could be any integer. Again consider the expansions of the superpotential W and R(ao) given 
by eqs.(3.13) and (3.14) respectively. On using eqs.(3.12), (3.13) and (7.1) in the shape invariance 
condition (3.9) and equating powers of ao one finds two sets of equations 

(i) n = pm, m = 0, 1, 2, . . . 

pm m 

9 pm (x) + ^9j(x)g pm ~j(x) = q m ^2g J (x)g m - j {x) - q m g' m {x) + R pm (7.2) 
j=o j=o 

(ii) n = pm + q, q = 1, 2 . . . (p — 1) 

pm+q 

9pm+q + J]] 9j( X )9pm+q-j(x) = Rpm+q- (7-3) 

These sets of equations are easily solved in succession to produce more solutions of eq.(3.9). 
Further, the energy eigenvalue spectrum can be easily obtained from eqs.(3.10), (3.14) and (7.1). For 
example, in the case where only Rx and R 2 are nonzero the spectrum can be shown to be (Eq = 0) 

= Jhj E^o/'- 1 ' + - Jby E^o) 2 ^ 1 ' , n = 1, 2, . . . (7.4) 

The energy eigenfunctions and the transmission coefficient for these reflectionless potentials can be 
written down immediately using eqs.(3.5), (3.27), (3.30) and (3.31). 

As an illustration, let us discuss the case p — 2 explicitly. The set of equations which follows from 
eqs.(7.2) and (7.3) is 

2m m 

92m( x ) +^9j(x)92 m -j{x) = q m '^2gj(x)g m - j (x) - q m g' m (x) + R 2m (7.5) 
j=o 3=0 

2m+l 

92m+iO)+ E 9j{x)g 2m +i-j(x) = R 2m +i (7.6) 

3=0 

and one can thus readily calculate all the g n (x). For example, in the case when only Rx and R 2 are 
nonzero it is easily shown that the first three g(x)'s are 

gx(x) = Rxx, g 2 (x) = (R 2 - qR x )x - ^Rfx 3 , 
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9z{x)^\R 1 {qR 1 -R 2 )x :i + ^-R\x\ (7.7) 
o 1 o 

Notice that we have chosen go(x) = so that again W(x) contains only odd powers of x, V±(x) are 
symmetric in x and SUSY is unbroken. The spectrum which follows from (7.4) (for p = 2) is Eq = 
and 

Ei^ = ^fog) 9 *" 1 + J , » = 1, 2, . . . (7.8) 

The g — + limits of the equations above again correspond to the one soliton solution of the Rosen-Morse 
potential, so that our results for a\ = qa^ can be regarded as multi-parameter deformations of this 
potential. Generalization to the case when an arbitrary number of Rj are nonzero is straightforward. 
Similarly, one can also consider shape invariancc in multi-steps along with the ansatz (7.1), thereby 
obtaining deformations of the multi-soliton solutions. 

Finally, consider solutions to the shape invariancc condition (3.9) for 

ai = — 7.9) 

1 +pa 

where < q, p < 1. We also assume that pao « 1 so that one can expand (1 +pa )~ 1 in powers of 
do- Further, assume that in eqs.(3.12) and (3.13) 

R(a Q ) = Ria Q + R 2 a 2 (7.10) 

and go(x) = so that W is again an odd function of x. On using eqs.(3.12), (3.13), (7.9) and (7.10) 
in the shape invariance condition (3.9), expanding negative powers of (1 + pa ) in powers of ao and 
finally equating powers of ao, one again obtains a set of linear differential equations. For example, to 
order Oq, the shape invariance condition looks like 

a g[(x) + al{gl+g' 2 {x)) = -(^-\^ x ) + f^^)\ g f g >( x)) + Rl a + R^. (7 .11) 

VI + pa / VI + pa / 

Expanding the denominators in powers of a and equating terms of order a and yields equations 
for the functions gi(x) and g 2 (x) which give 



Rix 

5i W = Y - ^' 92(x) 



pqR 1 . 
R2 + (1 + 9) 



(1 - q) x 



3 



(l + o) 2 (l + a 2 ) 3 

The energy spectrum which follows from eqs.(3.11), (3.15) and (7.9) is Eq = and 



(7.12) 



Et^RiY, q n +^E T H , n=l,2,... (7.13) 

fe[l+P«o(^)] fell+W^)] 2 

As usual, arL d r_ can be found using eqs.(3.5), (3.27), (3.30) and (3.31). Generalization to the 
case when arbitrary numbers of the Rj are nonzero is straightforward. Similarly, one can also consider 
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shape invariance in multi-steps along with the ansatz (7.9) and obtain deformations of the multi-soliton 
Rosen-Morse solutions. 

VIII. Summary and Open Problems 

Until now, the only known shape invariant potentials were such that the parameters a\ and do 
which appear in shape invariance condition (1.1) were related by a translation. In this paper, we 
have discovered a wider class of new shape invariant potentials for which a\ and ao are related by 
scaling, as well as in a variety of other ways. All these new potentials are reflectionless and have an 
infinite number of bound states. They can be considered to be q-deformations of the multi-soliton 
solutions corresponding to the Rosen-Morse potential. We were able to obtain the energy eigenvalues, 
eigenfunctions and transmission coefficients for these potentials algebraically. It was also possible to 
obtain analytical answers for the moments of these potentials, which should be useful since these 
potentials could not be explicitly expressed in a closed analytic form. The recently discovered self- 
similar potentials of Shabat and Spiridonov were shown to be a very special case of our shape invariant 
potentials. We were also able to obtain q-deformations of the one dimensional harmonic oscillator 
potential. This work has raised several questions which need to be looked into. Some of these are: 

(i) Just as we have obtained q-deformations of the reflectionless Rosen-Morse and harmonic oscillator 
potentials, can one also obtain deformations of the other simple shape invariant potentials? In 
particular, can one obtain deformations of potentials which are not reflectionless, say non-solitonic 
Rosen-Morse potentials of the form V(x) = — A(A + l)sech 2 a; for non-integer values of A. 

(ii) What are the various potentials satisfying the shape invariance condition (3.9)?. In this paper, we 
have significantly expanded that list but it is clear that the possibilities are far from exhausted. 
In fact it appears that there are an unusually large number of shape invariant potentials, for 
all of which the whole spectrum can be obtained algebraically. How does one classify all these 
potentials? Would such a classification exhaust all the known exactly solvable ones dicussed by 
Natanzon [19] ? 

(iii) The shape invariant potentials have been treated algebraically in this paper. An obvious interesting 
question is if one can also solve the Schrodingcr equation for these potentials directly? This should 
be possible, at least in principle. In that case the next question is if the Schrodinger equation gets 
essentially reduced to a hypergeometric or confluent hypergeometric equation or not. If not, then 
one would have generalized the concept of solvable potentials as introduced by Natanzon [19]. 

(iv) Now that a host of new shape invariant potentials have been discovered, it is worth asking if 
all the known exactly solvable potentials of Natanzon can be cast in a shape invariant form. In 
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fact one can ask an even more general question: can any exactly solvable problem in quantum 
mechanics (i.e. for which the Schrodinger equation need not necessarily reduce to a hypergeometric 
or confluent hypergeometric equation) be cast in a shape invariant form? In other words, is shape 
invariancc not only sufficient but even necessary for exact solvability, as first conjectured by 
Gendenshtein [6]? 

We hope to answer some of these questions in the future. 
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Figure Captions 

Figure 1 Self-similar superpotentials W(x) for various values of the deformation parameter p. The Taylor 
series are summed for x < 0.4 and the functions extrapolated from there by numerically solving 
the self-similarity condition, eq.(2.9). 

Figure 2 Self-similar potentials V-(x) (symmetric about x = 0) corresponding to the superpotentials 
graphed in Figure 1. 

Figure 3 A double well potential V-{x) (solid line) and its shape invariant, single well supersymmetric 
partner V + (x) (dotted line). The exact spectra are also displayed. Parameter values are n = 1, 
r 2 = -1, q = 0.3 and a = 0.75. 

Figure 4 Selected members of the one parameter family of isospectral potentials which includes the self- 
similar potential withp = 0.5. Note that a different choice of C\ (i.e. n) has been made compared 
to Figures 1 and 2. 
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